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ABSTRACT 


Peele vworetne f£luctuating velocity field in the 
developing region of a pulsating pipe flow is presented. 

Theoretically, three separate determinations of radial 
and axial perturbation velocity profiles were calculated 
involving varying degrees of analytical simplification. 

Heoctimentally, pulsating disturbances were created 
by a sinusoidal pressure wave exciter ins cele at the out- 
Pi—enea tes tack Circular wind tunnel facility. Cross- 
wire anemometer measurements of the pulsating velocity 
components compared favorably with the parallel theoretical 
results. Likewise, experimental determinations of the 
phase angles between velocity and pressure signals pone 


favorably with analytical values. 
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NOMENCLATURE 


Constants, coefficients of the polynomial 
approximation of the velocity profile 


Inside diameter of the pipe 
Paes lage 

Radial. coordinate 

Radius of the pipe 

Reynolds number 

Time 

Mean velocity in the pipe 
Maximum velocity in the pipe 
Axial component of velocity 
Radial component of velocity 
Pte teeCoorndinace 

Kinematic viscosity 

Density 


Angle of phase shift 
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I. INTRODUCTION 


The development of flow in the inlet region of a circular 
cylindrica] tube has been a subject of interest since the 
time of Hagen (1839) and Poiseuille (1841). Osborne Reynolds 
(1883) conclucted extensive investigations into the stability 
of pipe flow, and this area has become a classic problem 
with a sizable amount of literature published. However, the 
Senay Of fluid flow through a rigid tube under the influence 
of a periodic pressure difference is more limited. 

The theoretical study of pulsating, laminar pipe flow 
mimerimariay given by Th. Sexl [1] and S. Uchida [2]. 
Experimental data are limited to the works of Atabek, Chanc, 
and Fingerson [3] and the results of Florio and Nena [4]. 
With the introduction of the modern computer, several attempts 
have been made to solve the problem by numerical methods. 

The most recent work in this regard is that of Lew and Fung 
[5]. In the theoretical analysis the main problem has been 
the non-linear convective acceleration term in the Navier- 
Stokes equation. Uchida [2] was able to apply linear theory 
for the fully developed laminar pulsating flow. Lew and 
Fung {5] found this linearization to be valid in the entry 
region for high Reynolds numbers. However, this assumption 
leads to significant errors for small Reynolds numbers (the 


area of most interest). 





The present investigation was undertaken to study the 
fluctuating velocity field which arises in the developing 
region of a rigid walled conduit due to an imposed sinusoidal 
pressure gradient. Additionally, results were sought con- 
cerning any phase differences between the velocity and 
pressure Signals with a view toward later stability 


@ererminations. 





Il. THEORETICAL ANALYSIS 


A. SOLUTION # 1 OF PULSATING FLOW EQUATIONS 

For oscillatory flow in a circular pipe, the Navier- 
Stokes equations must be solved for u and v, the fluctuating 
velocity components, as functions of three variables (r,x,t). 
The velocities, wu and v, are independent of ¢, the azimuthal 


coordinate. The applicable equations from Schlichting [6] 


are: 
eee py et aL OF yy ou , Lau, aul (1) 
ot ox ox 0 Ox yr? = ax? | 
Iv Iv IV iL BE a7y lov, d’v vy 
Pree ens 5 jx \..2 °° ...° .2 22/'?? 

ox 1g rg OX r 

mee CONtinulty equation is: 
em OY a (3) 
OX oie ie 


Pic miosumcauthat uw = u + u, p = p + p, and v = v where 
u and Dp are the steady components and U, v, and p are the 
pulsating components. 

After linearization, two additional assumptions are made, 
namely, (1) that 6 is applied across the full cross-section 
Or eetewtibeyel.e., Dp 1S independent of r, and (2) that the 
gradient of Vv in the x direction is mild compared with the 


gradient in the r direction. 
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The resulting equations are: 


eu Vay ivy te (4) 
ox 12 Ghee tg 

sa | = Ot . ou __ 1 2 

aa Uae ve 5 ae (5) 

ou |, av, VL 

ax * or * x ° —) 


Suescertucion Of (6) into (5) yields 
Ju : 
TV 9x (7) 


The left side of equation (7) is recognized as that treated 
by Uchida [2] and the right side of (7) is known if ¥ and u 
ee enOwnN - 


Pie) 


By assuming v = R(r)X(x)e equation (4) may be 


solved by separation of variables. This technique yields 





toe - Se - = - w (8) 


Sil < 
w 
_ 
No 
<S 


Me Ce (9) 


and the R solution after assuming a series representation 


becomes 
pee a Br* | 8a-B a i B {20 _ 22 6 
“ae 9 192 192 |3 r 
, 2 [o% _ 082, 8° | 8, ie 
192 10 8 8 ee 





Therefore the complete solution for v is 








le 2 ay, 
v= Car eee If i a a oS x4 
| (11) 
2 D: 3 
_B {| 24 _ 2 Cee oO. os Cin eee 
"Dee ve B r* 92 |i07 48s * go/7 + ” 
where 
-0*U_ 
Qa = z (12) 
VI 
WU -iN 
oS a (13) 


and iG 1s the constant introduced in the separation of 
variable technique. Now that a solution for v is known, 
equation (7) may be solved. By substituting equation (11) 


gmto (7) and further assuming 


~ ? 
= A 
0 == Jo ae bir + bo 16 (14) 


1g nip etie 


_~L Xp. nN. 
; ae v Pne pier 1) (7) e (alg5)) 


equation (7) becomes a zero order non-homogeneous Bessel 


equation of the form 





ae y= p+ ps eee (16) 
ic v dus r 

The relationship in equation (14) allows determinations in 

the developing region to be made. The relationships in 

equation (15) concerns an oscillating axial pressure gradient 


and a resultant oscillatory axial velocity fluctuation. 


ik 





The complimentary solution of equation (16) is 
U, = C.J (2) + CY, (2) La) 


where z = v-i Xr and A - /% : 
The particular solution is obtained by the method of 
variation of parameters (Wylie[7]). 


U = u 


p 1¥, FB 


2% 2 (18) 


where an Jy (2) and ve = YQ (2), and 





; “Y9 
4 = ——_-————. R(r,x,t) (19) 
pee 59 4. 
er : 
u.7 = ————— R(r,x,t) (20) 
Y¥i%¥2 “Y2¥) 
with 
-int - S —~ 
Mita, xt) = Py aS V — u wv — (2a) 


The denominator of uy” and u,” fomthe wronskian oF (yy 1Y5) 


(Abramowitz and Stegun [8]) 





2 
eee as 22 
Wly,/Y>) = (22) 
iiererore 
y il 
= 23 
ec _ AZ 24 
Us 5 Jy (2) Rear, x; £) (24) 
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Entering the formulations for u and ¥ into the expression, Lor 
R(r,x,t) and then substituting the resultant expression 


equations (23) and (24) yields 


; TP 
u, = Ee . os fz YQ (2) (25) 
2 

+ hy fz Yo (z) + ho [2 YQ (2) 

- 2d, Py, Y(z) - A, [24 Yq (2) 

- de [2° Yo CZ) ao [2° YQ (2) 

- ha 2” x, Yo (2). = ho fe” Yq (2) 

= do [2° Yo (z) + | 


wiere the Various constants are listed in Table I. 


Sc ounm dls he hag 


ie 
yy a et? [Hu , Jp (2) - (26) 
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meer Integration to obtain the expressions for Uy and Us and 
Substitution back into equation (18), the complete expression 


for UW is obtained. 


B= CJ, (z) + P eV 6 (27) 


N a T 


x 1p | 22 _ 
A,e E Jy (2) | YQ (2) 


+ ¥,(z) | Jy(z)] +2 abe 


3 

2 eee een |= (e957) 

+ 9 | Jy (2) = | YQ (2) 

+ rpe™® [2 (24-627-16) 

+ vger? < eS-ase%ea25n) + 225 | YQ (2) ~ | JQ (2) 
+ ie 2 (26-3624+2082"7-96) 

2 ager? £ gl har oree) 

> OASz) 8 any wae [aq (2 z [yQ(z) 
4 Agen? 2 (a8 sasfaas04c4-20, 064224147, 456) 


The constant C, is made equal to zero by the boundary 


3 


condition that u must be everywhere finite. The constant 


C. is solved for by the boundary condition that u = 0 at 


rn 1 


iS 





The difficulty involved with this expression is obvious 
Thus a second attempt at a theoretical solution was made to 


obtain a more manageable solution for the ease of calculation 
B. SOLUTION # 2 OF PULSATING FLOW EQUATIONS 


In this second phase of the theoretical analysis, 


assumption that v is independent of x is made. 
equation (4) becomes 


the 
Thus, 

35 322 See 

Mane 9 ~ 3 

Gite or 12 Ge i 


(28) 
Again, an oscillatory form for v (now independent of x) 
is assumed, namely v 


R(x) oi (Ntt+¢) 
Perna lation into equation (28) 


SubStitutwon sol stnas 
equation 


yields a first order Bessel 


tir -(-5+4 +2] R=0 (29) 
ng 18 v 
with the solution 
2) SW ees (30) 
Incorporating this new expression for v alters the results 
FOr u, and u, Seeetewparelcular Solytion of u, 
expression for us becomes 


The new 
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TP..e e 
A Ty 58 3 
uy roe 2 J, (2) — 7 2 J (2) J, (2) C31) 
4 2 2 
. r r Z Z r Z, 
eM Aue 232 = - 2 2 
va 5 3 5 59 3 Jq (2) 
2 4 2 
: ) ney r Z nN 2. 
ci) || eg 7 a me Qos 2 
te 5 hoo2 3 ho32 D J, (2) 


As J.) (x) and Y, (x) ere both cylinder functions, the 


expression for u. can similarly be expressed by substituting 


Jl 


Ys eve J. The constants oe are tabulated in Table II. The 


final expression for u thus becomes 


ie C J (2) + u 


5 Jy (2) + U, YQ (2) 2) 


il 
The immensity of u has been somewhat reduced. However, 
it was felt that additional assumptions would yield a very 


workable solution. 


C. SOLUTION ¢ 3 OF PULSATING FLOW EQUATIONS 

After obtaining the two previous solutions for U, a more 
Simplified expression was still desired. An order of magni- 
tude analysis for u near the wall and for u near the center~ 
line was made, these being the two main areas of interest 
during the investigation. These order of magnitude determin- 
ations did lead to the finding of a simplified solution 
applicable across the entire cross section of the pipe. 


The expression for v is 


= Cc ai (Ntt?) 


1 J, (2) (33) 
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Pas eomonmiENeos boundary Conditaon that v = 0 at r = r 


i. (Nt+) w37 2 = | 
R Cy e J, (1 AX g) = 19) (34) 
which is satisfied either by Cc, = 0 or by 
ae?) i Gs Arp) = 0 (35) 


Ne@decting the trivial solution, then 


qi (Ntt+?>) 


R ea? Arp) = 0 (36) 


which leads to 


cos (Nt+¢o) ber, Ar, - sin(Nt+o) bei, Ar = 0 (37) 


ul 0 


or Simplified 


ber, (Xp) 
tan (Nt+¢) = bei, (x,) (38) 


where XQ = AX g- This form suggests a phase angle dependence 
upon both time and frequency. For t=0, the dependence of 
phase angle upon frequency of oscillation is: shown in Figure 
ie the constant Cy in this eigenvalue formulation is 
indeterminate because of the linearity of the problem. 
However, aS in all such cases, it is the form rather than the 
Mipltedaew or ene perturbation velocity which is of interest. 
With ¢ thus determined from these eigenvalue considerations, 


ee ai (Ntt?>) 


v is of the form V J, (2) where the arbitrary 


constant has been set equal to unity for calculation purposes. 
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Determinations upon u arise from equation (7) 


en ne ey (39) 


fe-ethe wall, the velocities u and ¥ appearing on the right 


Side of equation (39) approach zero. The velocity gradients, 
gu av a. 
ar and ie ir EIOuGmielarge, are Tinite. Thus, near the wall, 


equation (7) becomes that treated by Uchida, namely, 








2 ~ Ihe be ao 
iy u D OX at = 0 (40) 
Mien resultant solution of (40) with x = Ar is 
Meee a>’ x) 
a = oho 2208 (ee (41) 
WALL ig a eae 
0 0 
Near the centerline, u is maximum, V is zero and = 
cel ANTI 
approaches a value of 5 . Likewise, here ra equals zero 
3% cer 
and ay approaches a value of 5 . Thus, near the 
centerline, equation (7) reduces to 
2. _19P _du_ _g 2 2 
va po ox ot ~ Umax “1° oe 


The solution of this non-homogeneous Bessel equation is 


. ids 
x) Craie~  U 
eae ali MAX (43) 


0) PN 


rg 





limecies expression for u, the last term is negligible 
Compared to the first two. Thus, a complete solution for Ut 


applicable to the entire pipe cross-section is 
Le — ae — 1 (44) 


Ges last form is equivalent to Uchida's final result for 
mepeowaxtal fluctuation velocity component. Additionally, the 
companion normalized (Cy = 1) radial fluctuation velocity 


component is 


- ai (nt+¢) Jl x) (45) 


il 
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IIL. EXPERIMENTAL ANALYSIS 


A. EXPERIMENTAL APPARATUS 

The experimental trials were conducted in a horizontal 
Lucite tube 1.5 inches I.D. and 45 feet (300 diameters) long. 
A centrifugal fan was used to pull air from the stagnation 
conditions in the plenum chamber through the pipe and 
pulsator assembly. The inlet chamber consisted of a rectang- 
ular box three feet square in cross-section and four feet in 
Wenger, the inlet chamber was fitted internally with five 
6-24 mesh wire screens to eliminate inlet turbulence effects. 
Joining the plenum chamber and the tube was a two foot long 
settling chamber. The well rounded inlet horn was a quarter- 
Seetrom of a 9 inch by 18 inch ellipse molded of fiber glass 
and polished smooth. The tube, rectangular inlet chamber, 
and settling chamber were leveled and aligned to eliminate 
Gay eecurvarture effect. 

The pulsator assembly consisted of a sinusoidal oscillating 
proton arrangement driven by a MB 1250 Exciter. The pulsator 
was located at the end of the Lucite pipe immediately prior 
to the centrifugal suction fan. The piston arrangement 
consisted of a 12 inch diameter disc reciprocating over an 
11 inch diameter opening to the atmosphere within the 13 inch 
cubical outlet chamber. The pulsating pressure signal was 
created by vibrating the piston in a direction perpendicular 


to the direction of the main flow velocity. 


len 





Velocity measurements were performed with a 1248-L-T 
1.5-6 Thermo-Systems Inc. type cross-wire multi-sensor probe 
pecdeassociated electronics. The hot wire probe could be 
miserted at measuring ports located at 16, 20, 32, 40, 48, 
56, and 272 diameters downstream of the entrance. 

In all cases, the hot wire probes were mounted in a 
transverse micrometer mechanism, which was capable of placing 
ee escnsor to within .001 inch of any desired radial position 
except at the wall. Probe configuration limited placement 


meee wall) to within .050 inch from the rigid wall. 


Pee BAPERIMENTAL PROCEDURE 

The developing mean velocity profile was monitored at 
Eeacting ports in the Lucite pipe located at 20, 40, 56, and 
272 diameters downstream of the entrance. The quantitative 
data was recorded at these four stations for Reynolds Number 
values (R = Uuaxt/Y) of 2848 and 5398 respectively. 

The fluctuating pressure differentials for high and 
medium mass flows were obtained using BRUEL & KJAER Z inch 
condenser microphones (Model 4136). The mean pressure 
differentials were obtained uSing micro manometers. 

To measure the pulsating velocity components, the cross 
wire probe was installed at station 1, twenty diameters 
downstream of the entrance. Radial distribution of velocity 
amplitudes and phase shifts were taken at Sinusoidal pressure 
oscillations of 15 cps and 20 cps and at high and medium 


mass flow rates. All time-dependent data were recorded on a 
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two-channel Sanborn Model 296 recorder. The physical setup 
Se the electronic equipment utilized to measure the OSei1 Tiating 
velocity components and pressures is shown in Figure 2. 

On the two-channel recorder v and U data was recorded 
alternatively; the input pressure signal was continuously 


Memicored. 


feeeeoaAreRIMENTAL RESULTS 

Initially, mean flow velocity data was evaluated at station 
1 where the oscillating velocity determinations were to be 
obtained. A smooth curve was fitted to <he mean flow data by 
the least squares method. Figures 3 and 4 illustrate the 
computed pilots for flow Reynolds Numbers of 2848 and 5348 
respectively. Table III gives the numerical results for the 
developing velocity profile at the succeeding stations. 

Secondly, verification of the phase shift, ¢ , inherent in 
the radial velocity expresSion was sought. After choosing 
an arbitrary t = 0 setting upon the records of the pressure 
and velocity Signals, a subsequent marking of parallel times 
allowed for a direct verification of analytical values as 
shown by the representative tracings in Figures 13 and 14 
for a pulsation frequency of 20 cps. While equipment limita- 
tions forbid a conclusive verification at this frequency; 
nevertheless, the appearance of the 90 degree phase shift 
between the axial velocity and pressure signals on the center- 
line.which is evidenced in Figure 15 and which is also 


predicted by theory makes the radial results less tentative. 


758) 





Finally, a comparison of theoretical and experimental 
meeues for both radial and axial fluctuation velocity 
components VaomceteniitlieC ds "“oteowmein Pigures 5-12. In the 
case of the axial velocity component, the comparison between 
experimentally measured values and theoretically calculated 
values is direct whereas for the radial velocity component 
the theoretical curves incorporated an arbitrary amplitude 


factor. 
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IV. CONCLUSION 


The final theoretical conclusions follow from the anamolous 
physical situation represented by equation (28) which demon- 
strates the viscous decay of the radial velocity component. 


Removing the x dependence of v, removed the production term 


-—- 9 Fe Se ee 
(u = which was necessary to generate the v term initially. 


Thus, subsequent theoretical determinations are made upon a 
velocity component which is not strictly allowed by the linear 
analysis pursued. This anamoly is represented by the | 
appearance of the arbitrary constant in the final expression 
for the radial velocity component. 

A second manifestation of this anamoly is evident in the 
independence of time that exists in the final expression for 
V since 


ee Sil 
(nt+¢?) = tan ber, X_/bei, x 


0° 
fimcmtehe initial form of v which follows from this result must 
persist. Both the form and persistence (since RMS evaluations 
are involved) seem to be suggested by the comparison of the 
theoretical and experimental values for the radial velocity 
component portrayed in Figures 9-12. Again it is to be noted 
that the theoretical curves here incorporate an arbitrary 


amplitude factor. However, as in all linear eigenvalue 


IES 





problems, it is the form rather than the absolute value of 
fie Cigenfunction which is important and which is verified 
experimentally here. 

The comparison of values in the case of the axial velocity 
component as expressed in Figures 5-8 iS more conclusive. 

This data, though limited in nature verifies previous (Richard- 
son [9]), theoretical and experimental determinations and 

thus helps substantiate the parallel radial data which was 
Simultaneously recorded. 

These mean determinations portray another characteristic 
of pulsatory flows, namely, that the maximum mean flow does 
not occur at the center of the pipe but in a region near the 
wall. The location of this region approaches the wall as 
pulsation frequency increases. As evidenced in both the 
radial and axial data, ce high value of the argument in the 
Bessel functions involved in the velocity expressions makes 
the wall area values predominant. Thus the central portion 
of the flow is recognized as being practically unaffected by 
pulsation. Moreover, the value of axial fluctuation velocity 
iS a non-zero constant on the pipe axis : another result 
wetmrred by the data in Figure 5-8. In contrast, the zero 
value of radial velocity fluctuation on the pipe axis is 
evidenced by Figure 16 which portrays the radial velocity 
component recorded with the same gain readings as pertain to 
all other tracings shown. The appearance of this figure is 
that of a filtered noise signal and demonstrates the favorable 


signal to noise ratio evidenced by the other recorded data. 
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In addition to the results concerning the existence of 
a radial fluctuation velocity component, the second new 
result warranted by the accompanying analysis was the appear- 
ance of a frequency dependent phase shift in the formulation 
of this velocity component. Results concerning this phase 
shift are shown in Figures 13 and 14 where determinations 
for t = 0 values approach those values of phase shift 
suggested by the analysis for a pulsation frequency for 20 cps. 
Equipment limitations did not allow verification of non-zero 
time values. However, as stated previously, the appearance 
Of a 90 degree phase shift in the axial velocity and pressure 
data recorded in Figure 15 tends to substantiate the companion 
radial velocity data. 

In conclusion, the limited experimental program reported 
here only demonstrates trends toward the parallel theoretical 
determinations. A more comprehensive experimental program 
Peercderequire both new electronic filters with band pass 
capabilities in the .0l1 cps range and a new pressure monitoring 
device to replace the microphone system which has severely 
attenuated response below 20 cps. In view of these require- 
ments, a re-evaluation of analytical results is recommended 


before a more extensive experimental phase is commenced. 
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Aes dois a 


Constants for series solution of v; solution # 1 of the 
Pulsating Flow equations. 


—W xX 
7 Ci ayboe T 
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Vv 
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Cjayte 
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Pe LE iT 


Constants for v independent of x solution; solution # 2 
of the Pulsating Flow equations. 
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1. hoy eS) oe 
Lo Ynv 
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TABLE IIf 


Coefficients of the polynomial approximations to the velocity 
profiles of the form: 


Mees 4g ptt \e pb. [© ar ai 
Um 0 1 Lo 2 Lo 3 Lo 
x/d mass flow rate Gr ae ED 2 
20 high Oo ieee 0964 I SSS 3: 2 Oe 
medium O27 G0 oa o ele oe SNe Cay aS 
40 high MOS = 2G 2.35016 ens) He 
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Figure 1. Phase Shift (degrees) vs Frequency (cps) of 
v compared with pressure signal at t = 0. 
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Developing Velocity Profile; x/d=20, 
R=2848. 


34 














r/X 9 


Figure 4. Developing Velocity Profile; x/d=20, 
R=5348. 
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Pagume 5. Axial Velocity Component; R=5348, 


Peeq = 20 cps. 
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Figure 6. 


Axial Velocity Component; 
Freq = 15 cps 
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Figure 7. Axial Velocity Component; R=2848, 
Freq = 15 cps. 
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Axial Velocity Component; R=2848, 
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Figure 9. Radial Velocity Component; R=5348, 
baege= 20) cos. 
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PegicestO. kaaial Velocity Component; R=5348, 
Pwear = 15-Ccps. 
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Figure 12. Radial Velocity Component; R=2848, 
Heede= 15° cps. 


43 








a+ a = { eee * 
pee fae . =, 
we oes + —e—— + Atay py = 
. \ a 
| ee i 


an 
en aa RR pi Se gS mi 
a - 9 S'* Fed & = 

A e. 


eee (co 
pa 


eae 









ee 


i 
| 
| 
i 
t 
| 
: 
! 


. 
eee | 


~~ = 
+¢ + 
oc ow 


' 





ero At ith doce ene 
Pe comes sy . "™ 
—— ; agra 
+ 
ge SL : rh ge 4 . . o...,, bs 
i + crm ewe ce 2 ee ee Sem + 6 ae ora = & Vem 


i —.* + 
7) ’ ‘ e- 


een ean 8 


| 
a 


: | \ 
iV. 
| 


a 
e 
La 
= he 
° 


y 


; ; j eee a 
‘ ' ¢ 4 , Alas ' fess 
_ 4 a a 5 


a ir ony 


J 
a 
] 
| 
‘ 


\ 
‘ 
1 
‘ 
\ 
t 
rs 


| 
a 
I 
tot 
i" 
| 

| 

| 
Bz 
7 
; 
Leo 
4g 
[3 
| 
oF 


4 
nf ¢ j 
eo: eee ee eee we —e o> ~6 - ee emer bem ° a e-eatiatbe - SO a ee 
nr er te at et ee ' i : 
' oo eo te ae 4 P weenwe i 
\° . t era o- of . e ‘ ‘ 
Ree A Eee pee 0 ‘ rear ; 
aa ene ee wd At 28g . 2+ t= ee | oe cee * bees @ we eve 
erceper . ese ‘<n + © meee tg ai ‘ S a 
owe oe 6 es 8 5 Sa | { pve eens ’ ° 
seve f be wo heecs+. ee Fe ’ : (ws ’ “ 
bo be pee o whew ent. wemerewe -' t - : . aan ° a Sere yee = 
SS a ee aes eee te 0 eo eet me mee 0 em my 9 oo + eee eee Ot -  -e ae -2e- 
mame yb T a! 
Ome ry, 1 4 ‘ ' ae ee es 
o a ' 4 1 ‘ eae - t +e -G-eret 
al ’ . ' ‘ A - ’ t ++ Me ene 
eos . : r- 1 t <2 1 ' to e+e 
= ee — a ee rr ne tee ah ee eh ss apie | = eo 8 eee 
en nf . ! ‘ . ' a . a = 
. { . + = . . t °¢ 
oie Sr = 4 > * eee. ‘ . . eos 
a Fr es Haw as ae = ° -- ame § eeu —@ t= -¢ onaeem -e e+ 6 s . -et oe = « &eress~ 
ee a - . tee ’ + . 
—————— sewed ie ek 
. 
4 7 pe 1 ' «.., ‘ Pa + . 
Fac . > . © tuhcd ae ‘ e-« 
ww > bes. — - BRSo Sc - ~ we 9 ee Dae Oe gare gS ee | $ a 
S++ wee ewe. ’ ° o- 





Ft * + eee we fy + . 


: 
‘cee ete) eens ; ® S ‘ . ! ieee ‘ “APs, : is 
je ee nea a ae ee Beg nr eee ee) ae eS es : 

¢ e- . 

: : I 





; . 
’ ae ate 4 — + be ry t- Soe ey, 
ee eee tcc - ~~ . == — =) -eere-bee awe . -- + o wb - . 

. . ° . - py enema | . -o* ’ ° - 
’ AGS See | a 2 < ren . Seca 
i - - Nees . { A ae . 75 Qn ~ eee oe 
' -* ot. 


- oe oe | - - -* ° -~ $+ ot wee 7, ry 
re ee a ee ee $$$ rrp tre 2 de ete ee ee ee SS Se ee ee 





° Hae | . eae r 9 8 tis 1 - 4 - + . ‘4 reer . e-e2 ee 
meme feocee . ’ . ' e —e 41 ee 
Mipieieh te a ——_—” : ‘ (ake : b : SE 
t- -$ ot ’ s . 2 + 0 om é cc 
ee ee ee gy eee ees Se ee OCs eae conical 
. - ece ° 4 mat) OF tees, ' . A t--o~ = ars e eee 
- Sa] ee ’ See . ree) j Poaceae Ano 
Eee Sea ee : me a ee. Ss fe ! ; ee 
3 ate BMt at =. . at oe l 3 chs At be ec og ‘ a i 
-—— 2 a om- = - = a i i ee ee --—= beng —-+-—- ee ee 
ae Tine hel, Ss ne ’ is 8 ee pe, i * . oS re —ooe qe o< ~~ + 
neo) Ve eS 4 a \ Het sf : 5 yn 
mes - bee eee a j . t [ <tes ‘ . + +e +o 
Rea meres e ‘ ean t : t We eerie “hes eo gee ae as 
Ce a a a ee ow - cre- mm at mp ee 2 ee es oe ee see a a re a dy 
_- by ea) caine - ‘ a ; t . * —~— ‘ <a - 4 
eae : - ; ss i : . Re oo ; : yO Lea. 
' e@rrcesce ’ + . ’ Co ot elon i a << . ° + 0 om 
eorere be * . ' 7? . - e ry ams -; a . - ee | 
ee -——— 6+ os, -- ee ne ee ee re ge ow eee - ee Se | a on 
e- . - ‘ ' . ‘ ws 0g . ’ j . ' ~— — -~e¢ o- 
SS Ue a} ais e $ aa ' * co ' — a6 4 Baie 
ie <'» a ar a ao FA aa es Fa ee) eee) ra ss 
ip | I yee) Gar jk —_— me = 06 wom. we — Oe pars Pa a aaa necpee ea a 26 he ee he ee mh a ee 


Figure 13. Recorder Tracing of p (left) and 


v (right) Radius = 0.7, R=2848, Freg = 
POmGos ye @ < 45° 5 


44 





eS SCs eee ee? Se 

oe oe ee es 

ASH lea es Cepia 
es - 


ee ee 







ee pe emt eee 2 ay 


a =f: 








aan A er @ wae coe SE) amen (ee a + > © ale oN) = «= 
eee ee ee ie cee ee ee ee eee oo me oss coeetoon 

= 
= anes + . . 


Set ee ee ar eae ee oe 
a eae en amhamiadl 


eae ee Peer Seer 


i es 

















ee a, 


s-ae 














ie) toes @ wm er d= 6 cc eh eee ge == —- © 


Tree ceo 








- - ne a 






. 
eee we oe 0 0-6 oe cb oh 8 eee = Py 
* 


op eae es 
Ce ee 






re oa ae eee Tesrapeee eee | 


r: 





+ wee Powe 
eee 


4 


Se ee ee ee 


- 
= 9° , 
= ‘ 
. . 

- ct ee 















- @ere tee rune Ce ee oe 
. : ’ +4 =< a6) 6 6 elite + > ares e ok we 
OfP OT) 0 DOOR RIE Or ss a ae Se Det = — +44 in be Se eee apie = er ao megs 
ttn be eee - oe = ¢ e+ eter ge ~ *eeege . eo + 7+? e+ tee bee 
cs epee UE — oe ee be ae —— 
A eae - delaras oe - | tf wrvce ee 
° oe Gewese ee fry . + cee - 
° ' +eeeote - oa Wee t Pegi. tae be 
++ ee+-te + 4 cote te a ce 7 

tee 






-{ 
rs 





- 2 0-4 terete 










ee ee ed 
+h wc ees rove ese 


Ra aA ee 
me ae moro aras 
















Se eee oe ee 
Cl bony eat Sree Cy 
a ee 2 





e--- cp eces 
¢ e+ e¢eewte 

‘ SK oa + iw 
' wmoee- 






a Si 2 oe ee 
women eet tet be ree bee - 


| - -- @e = 0 we rteel + 








c= ere “ee tew es fy o ’ 
° Ls a er+pees-e { 
ee ee ee ee ee - tee Roel | 
ee ee ree — 
= Seo) 4 eee As le Os 


Ae CRO Se Ae more ner cmomc Ie : 2 

Ceo opt Fa Oe eae (het ee ee ere 
Ce ee oe en or nas om Sit ee ¥ ee Ebay Gamer in 
So ee oe meee ene emg —- --- = bererer = eof) woke dt + ee wwe bam wes § - ote Hee spe 


a e+e pate 











eee 
a) 
oe 
4 
° 
° 
o- 
rh 
ben 
om 
“ 












. eevee -tinir ee ce oy 
Zea rasta Ge ene ny ect cee ed, 

a Ral a (ee ee peaee 
oe we b-o- . Sok, hello Cag Ciseo ben 
ea ee ae op pen weunn EETe a 
2 ee eee we ees 669 2S SS ae ee 7 ae a ee ane 
= =) eee ee we Se Sas SUI CEN pM Oe. OT, aie = 0 oho 60 ere-4 one = tr - Ge -ombmend 
m PL Go +3 a 


(SSS RS ae 







a ee 





Figure 14. Recorder Tracing of p (left) and 
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v (right) Radius = 0.7, R=5348, Freq = 
20 cps, ¢ * 45°, 
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Figure 16. Recorder Tracing p (left) and 
VeceraieleRadzus — 0.0, R=5348, Freq = 
lomecosemy = 0. 
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